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RADAR  BIAS  ERROR  REMOVAL  ALGORITHM 
FOR  A  MULTIPLE-SITE  SYSTEM 


1.0  INTRODUCTION 

Recent  investigations  as  to  the  feasibility  of  operating  a  multiple  platform  sensor  integration  sys¬ 
tem  [1-3]  have  led  to  the  problem  of  removing  bias  errors  from  transmitted  data.  Bias  errors  are  those 
errors,  inherent  to  each  sensor  system,  which  may  have  been  introduced  during  the  construction  or 
alignment  of  the  sensor  or  are  present  as  a  result  of  equipment  failures.  With  three-dimensional  radar 
systems  the  problem  reduces  to  one  of  detecting  the  bias  errors  that  may  be  present  in  the  range 
azimuth  and  elevation  measurements  made  by  each  of  the  radars  within  the  network  of  participating 
platforms.  The  solution  of  this  problem  is  fundamental  to  the  successful  operation  of  a  multiple  plat¬ 
form  sensor  integration  network.  If  bias  errors  cannot  be  detected  and  their  effect  reduced  to  a  toler¬ 
able  level,  the  advantages  to  be  gained  from  the  exchange  of  information  between  platforms  cannot  be 
realized.  As  a  first  step  in  this  direction  an  algorithm  has  been  developed  for  the  special  case  of  two 
fixed  radars  with  known  locations.  The  algorithm  has  been  tested  with  simulated  data  and  with  real 
measurements.  The  measurements  were  made  on  targets  of  opportunity  by  two  radars  located  at  NRL’s 
Chesapeake  Bay  Detachment  (CBD)  and  at  the  Applied  Physics  Laboratory  (APL)  of  the  Johns  Hop¬ 
kins  University.  The  results  are  presented  in  Sec.  3.0  and  4.0.  With  minor  adjustments  to  the  algo¬ 
rithm,  the  same  basic  concept  should  be  applicable  to  the  moving  platform  case. 


2.0  ALGORITHM  DEVELOPMENT 

The  problem  which  was  taken  under  consideration  is  described  as  follows:  Given  the  time  history 
of  measurements  (range,  azimuth,  and  elevation)  on  several  distributed  targets  as  measured  by  two 
separate  radars,  determine  the  measurement  biases  for  each  radar  under  the  assumption  that  the  exact 
location  (latitude,  longitude,  and  elevatK  '  of  each  radar  on  International  Geoid  is  known. 

A  two-dimensional  representation  of  the  geometry  of  the  problem  is  shown  in  Fig.  1.  Coordinate 
systems  are  established  at  each  site  with  the  x-axis  pointing  due  east,  the  y-axis  pointing  to  true  north, 
and  the  z-axis  (not  shown)  pointing  outward  along  a  line  joining  the  site  and  the  Earth’s  center. 
Azimuth  is  measured  clockwise  from  the  y-axis,  elevation  from  the  horizontal  plane  containing  the  x 
and  y  axes,  and  range  is  the  Euclidian  distance  between  the  site  and  the  target.  The  range,  azimuth, 
and  elevation  measurements  can  readily  be  transformed  to  the  local  rectangular  coordinate  systems  at 
each  site.  This  allows  us  to  express  the  measured  position  vector  of  a  target  with  respect  to.  Site  lias 

AT]  =  X\  +  A\B\  +  A\N\,  (1) 
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Fig.  1  —  Geometric  representation  of  the  problem 
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r  —  range,  a  —  azimuth 
’  and  7)  *=  elevation 


with  Xi  being  the  position  vector  of  the  target  at  some  instant,  Xx  being  the  mean  or  true  position  at 
the  same  time,  Bl  being  the  vector  of  bias  errors,  Ar1  being  zero  mean  noise  on  the  measurements,  and 
Ax  being  a  matrix  of  partial  derivatives  which  are  derived  in  the  appendix.  Equation  (1)  is  in  effect, a 
truncated  Taylor’s  series  expanded  about  the  target’s  true  position.  The  measurements  made  at  Site  2 
can  also  be  converted  to  the  local  rectangular  coordinate  system  at  Site  2  and  transformed  to  the  local 
rectangular  coordinate  system  at  Site  1.  This  enables  us  to  express  the  position  vector  of  the  target  as 
measured  by  Site  2  with  respect  to  Site  1  as 

Xi  —  +  A  2B2  +  A  2N2,  (2) 


where  A2  — 
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i.e.,  A 2  is  the  matrix  of  partial  derivatives  representing  the  rate  of  change  of  Site  l  rectangular  coordi¬ 
nates  due  to  changes  in  the  measurements  made  at  Site  2,  B2  is  the  vector  of  bias  errors  at  Site  2,  and 
N2  is  the  vector  with  elements  that  represent  the  noise  on  the  elements  made  at  Site  2.  If  the  measure¬ 
ments  at  Sites  1  and  2  are  not  made  at  the  same  instant,  it  is  possible  to  predict  a  "measured1  position 
at  some  common  time  and  still  have  Eqs.  (1)  and  (2)  hold  provided  the  partial  derivatives  arc  not 
changing  significantly.  When  Eqs.  (1)  and  (2)  represent  coincident  measurements,  X}  -  X2  and  it  is 
possible  to  subtract  Eq.  (2)  from  Eq.  (1)  to  yield 

A*=  AtBj-  A2B2  +  AlNl-  A2N2,  (3) 

where  lAA'I  represents  the  Euclidean  distance  between  the  two  measurements.  Equation  (3)  can  also 
be  expressed  in  the  form 
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where  M  =  \A  j. 


A*  =  MB  +  N, 


Bi 

Nx 

B  = 

b2 

,  and  N  —  M- 

n2 

(4) 


Equation  (4)  is  now  in  the  form  of  the  observation  equation  for  linear  estimation.  In  this  case  the  AX 
vector  represents  the  measurements  and  the  B  vector  of  bias  errors  is  the  state  vector.  The  state  equa¬ 
tion  is  given  by 


Bit  +  1)  =  /  •  Bit). 


(5) 


Note  that  the  state  transition  matrix  is  the  identity  matrix  since  the  bias  errors  are  assumed  to  '>e  con¬ 
stant  over  the  period  of  interest;  i.e.,  the  bias  errors  at  some  future  state  (/  +  1)  are  equal  to  the  bias 
errors  at  the  current  state  (/). 


With  this  formulation  (Eqs.  (4)  and  (5))  it  is  possible  to  predict  the  bias  errors  using  the  Kalman 
filter  algorithm.  (See  Ref.  4  for  a  typical  application.)  The  six  steps  involved  in  the  recursive  algorithm 
are  as  follows: 

Step  1.  Calculate  one-step  prediction, 

BCrlr  -  1)  =  BU-  1|/  —  t).  (6) 

where  the  circumflex  signifies  an  estimate,  B  (/!/  —  1)  indicates  the  predicted  bias,  and  Bit  -  Hr  — ,1) 
indicates  the  smoothed  bias  of  Bit). 

Step  2.  Calculate  the  predicted  covariance  matrix  P(t\t  -  1)  front  the  smoothed  covariance 
Pit  -  Ur  -  1)  bv 

Pit\t  -  1)  =  Pit  -  11/  -  1).  (7) 

(Note  that  for  this  case  the  predicted  covariance  matrix  and  state  vectors  are  assigned  the  values  deter¬ 
mined  at  the  previous  time  step  since  the  state  transition  matrix  is  the  identity  matrix.) 

Step  3.  Calculate  the  predicted  observation  AX (rl/  —  1)  by 

AXitU  -  1)  =  Mit)Bi'U  -  1),  (8) 

where  Mit )  is  the  measurement  matrix  Mat  the  rth  sample. 

Step  4.  Calculate  the  filter  gain  wit )  using 

wit)  «  l)M(t)  lMit)Pitlt  -  1  )Mit)  +  R  (/)]"',  (9) 

where  the  tilde  indicates  the  transpose  of  a  matrix  and  R  it)  is  the  covariance  matrix  of  the  noise. 

Step  5.  Calculate  a  new  smoothed  estimate 

BitU)  =  B{t\t  -  1)  +  mK/HAAT/)  -  AXitU  -  1)],  (10) 

where  AX  it)  is  the  measurement  vector  AX  in  Eq.  4  at  the  rth  sample. 

Step  6.  Calculate  a  new  covariance  matrix  , 

P(fl/)  =  [/  -  wit)  •  Mit))  Pitlt  -  1).  (11) 


The  measurement  covariance  matrix  R  is  the  covariance  of  N;  i.e.,  . 

R  =  cov  (AO  =  £{(/!,-  JV,  -  A2  -  N2)  (/l,  •  Nx  -  /f2,y2)r}.  (12) 
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Under  the  assumption  that  the  partial  derivatives  are  changing  slowly  <n  the  neighborhood  of  the  target 
location,  a  good  approximation  of  the  covariance  is 

cov  (AO  =  A\(N\Nl)Af  +  A2(N2N2)A2  (13) 

or  alternatively 

cov  (AO  =  M  (N  ■  NT)  MT.  (14) 

it  should  be  noted  that  in  calculating  the  partial  derivatives  which  are  the  elements  of  the  M  matrix, 
that  the  smoothed  estimate  of  the  bias  errors  is  applied  to  the  measured  range,  azimuth,  and  elevation 
to  provide  an  evaluation  closer  to  the  true  position  of  the  target.  The  matrix  ( N  •  Nr)  is  a  diagonalized 
matrix  with  the  diagonal  elements  being  equal  to  the  variances  in  the  measurements  of  each  radar. 
These  were  assigned  typical  values  for  the  radars  that  were  used. 

3.0  RESULTS  WITH  SIMULATED  DATA 

To  check  out  the  algorithm  described  in  Sec.  2.0,  a  set  of  simulated  tracking  data  was  developed. 
Stationary  targets  were  chosen  to  simplify  the  process.  When  moving  targets  are  being  tracked  further 
complications  arise.  Since  the  measurements  at  both  sites  are  not  made  simultaneously  a  prediction  of 
one  of  the  measured  positions,  to  some  common  time,  must  be  made.  This  prediction  process  is  hot 
required  for  stationary  targets  since  the  predicted  position  is  the  current  position. 

Six  stationary  points  in  space,  surrounding  the  locations  of  the  two  sites  (NRL’s  Chesapeake  Bay 
Detachment  and  APL)  were  selected  as  targets.  The  locations  of  the  targets  and  their  respective  lati¬ 
tudes  and  longitudes  are  shown  in  Fig.  2.  All  of  the  target  points  were  assumed  to  be  at  an  altitude  of 
10,000  m. 


N 
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Measurement  data  were  generated  by  selecting  samples  from  a  Gaussian  distribution  derived  fr< 
a  random  number  generator  and  by  adding  these  samples  to  respective  range,  azimuth,  and  elevati 
coordinates  at  each  site.  Bias  errors  that  would  adequately  test  the  algorithm  were  also  added  to  i 
measurements  at  this  point.  This  resulted  in  a  set  of  simulated  measurements  (range,  azimuth,  a 
elevation)  for  each  of  the  six  targets  as  measured  by  both  sites.  The  difference  in  the  location  of  ei 
target  as  measured  by  both  sites  was  used  as  the  input  to  the  Kalman  filter.  The  process  of  simulatin, 
detection  of  all  the  targets  by  each  site  was  repeated  50  times  or  equivalently  50  rotations  of  each  rad; 

The  computer-generated  bias  errors  as  predicted  by  the  algoritlim  are  presented  graphically 
Figs.  3  and  4.  They  show  the  estimated  bias  error  as  a  function  of  time;  the  time  units  being  i 
number  of  scans  completed  by  both  radars  from  the  initial  time  until  the  time  at  which  the  estimate 
made.  Several  cases  with  a  wide  range  of  induced  bias  errors  were  simulated.  The  case  which  is  illi 
trated  in  Figs.  3  and  4  had  the  following  induced  errors  at  the  APL  site: 

range  bias  error  =  2  n.mi 

azimuth  bias  error  =  — 3° 

elevation  bias  error  =  3 

and  at  the  NRL  site: 

range  bias  error  =  —2  n.mi 

azimuth  bias  error  =3° 

elevation  bias  error  =  —3  °. 

The  results  presented  in  Figs.  3  and  4  are  from  a  single  case,  i.e.,  both  sets  of  bias  errors  were  indu< 
simultaneously.  As  can  be  seen  from  the  figures,  a  reasonably  good  estimate  of  the  bias  errors  .can 
obtained  after  only  three  or  four  scans  of  the  radars.  The  curve  that  represents  the  estimate  of 
elevation  bias  is  not  as  smooth  as  the  other  curves.  This  is  due  to  the  value  assumed  for  the  variai 
in  the  measurement.  The  standard  deviations  assumed  for  the  range,  azimuth,  and  elevation  measu 
ments  were  400  m,  0.5°,  and  1.0°,  respectively. 
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From  the  preceding  results  it  appears  that  bias  errors  can  be  routinely  determined  when  several 
targets  are  dispersed  about  a  pair  of  widely  separated  (35  n.mi)  radars.  However  it  must  be  remem¬ 
bered  that  this  represents  an  idealized  situation  in  which  the  targets  are  stationary  and  evenly  distri¬ 
buted  (see  Fig.  2).  Because  of  this  consideration  and  the  availability  of  real  data,  collected  by  both  sites 
on  targets  of  opportunity,  further  testing  of  the  algorithm  was  carried  out. 

4.0  RESULTS  WITH  REAL  DATA 

During  September  1979,  the  SPS-39  radars  at  NRL  and  APL  were  used  to  simultaneously  record 
detection  data  on  targets  of  opportunity.  Subsequently,  these  data  were  processed  by  APL  to  identify 
tracks  and  develop  track  histories  from  each  radar’s  detections  and  to  correlate  tracks  of  common  tar¬ 
gets.  APL  also  utilized  their  tracking  algorithm  to  predict  the  position  of  all  tracks  at  the  time  of  north 
crossing  of  the  APL  radar;  i.e.,  regardless  of  when  detections  were  made  an  estimate  of  the  positions  of 
all  targets  was  made  corresponding  to  the  time  at  which  the  APL  radar  had  most  recently  swept  by  true 
north.  The  estimated  position  of  each  track  at  that  time  was  given  in  rectangular  Cartesian  coordinates 
centered  at  each  site  and  it  was  in  this  form  that  the  data  were  made  available  to  NRL.  Fortunately  the 
data  in  this  form  were  amenable  to  the  algorithm  described  in  Sec.  2:0.  All  that  was  required  to  make 
it  directly  applicable  was  the  transformation  of  the  NRL  tracks  from  the  NRL  coordinate  system  to  the 
APL  coordinate  system.  For  each  scan  of  the  APL  radar  each  correlated  track  was  examined  to  see  if 
both  radars  had  detected  the  track  since  the  previous  north  crossing.  If  this  proved  to  be  true  the 
difference  vector  (see  Eq.  (4))  of  the  predicted  positions  was  used  as  the  input  to  the  bias  removal 
algorithm.  For  the  data  that  were  used  this  could  mean  as  many  as  15  updates  by  the  Kalman  filter  in  a 
single  scan. 

The  effect  of  applying  the  algorithm  to  a  specific  track  can  be  seen  by  examining  Figs.  5  and  6. 
Figure  5  is  a  plot  of  the  tracks  developed  by  both  radars  for  a  common  target.  The  cracks  are  displayed 
in  the  horizontal  plane  with  APL  being  located  at  the  origin  and  NRL/GBD  located  approximately  35 
n.mi  to  the  soutneast  of  APL.  The  target  is  approaching  APL  from  the  east.  The  result  of  correcting 
the  t>acks  by  amounts  equal  to  the  estimated  bias  errors  is  shown  in  Fig.  6.  In  this  case  the  start  of  the 
track  was  coincident  with  the  start  of  the  bias  error  estimation  process;  consequently  initial  estimates  of 
the  bias  errors  were  applied  at  the  beginning  of  the  tracks.  The  estimated  values  from  both  radars  that 
were  used  to  correct  the  tracks  are  also  plotted  in  Fig.  6;  however,  using  the  initial  values  of  the ‘bias 
errors  to  correct  the  tracks  may  not  be  representative  of  a  working  system.  A  conceptual  system  would 
probably  have  a  fully  developed  set  of  bias  errors  available  for  correcting  the  tracks  or  more  likely  the 
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raw  data  and  updated  estimates  of  the  errors  would  be  made  on  a  periodic  basis.  To  assess  the  effects 
of  using  a  fully  developed  set  of  bias  errors  the  track  data  were  processed  for  15  min,  using  the  algo¬ 
rithm,  and  the  following  set  of  bias  errors  was  obtained: 

on  the  APL  radar  0.224  n.mi  in  range 
2.033°  in  azimuth 
0.119°  in  elevation 
on  the  NRL  ;adar  0.416  n.mi  in  range 
-0.691°  in  azimuth 
0.505°  in  elevation. 

These  errors  were  applied  to  the  tracks  of  the  target  shown  in  Fig.  5,  and  the  results  are  shown  in  Fig. 
7.  Again  the  results  are  quite  dramatic,  and  a  significant  improvement  over  Fig.  6  can  be  seen  in  the 
early  stages  of  the  track.  Undoubtedly  if  such  information  were  available  to  make  corrections,  it  would 
simplify  the  correlation  process  and  greatly  reduce  the  number  of  false  tracks  in  a  multisite  sensor 
correlation  system. 


Fig.  7  —  Tracks  corrected  with  fully  developed  bias  errors 


To  further  validate  the  algorithm  it  was  decided  to  apply  the  technique  to  a  target  in  another  qua¬ 
drant,  and  a  target  was  selected  from  the  southwest  quadrant.  The  target  started  at  a  pqint  approxi¬ 
mately  40  n.mi  west  and  60  n.mi  south  of  the  APL  site.  Initially  it  was  headed  toward  th«  NRL  site 
but  near  the  end  of  its  flight  it  veered  toward  the  APL  site.  The  two  tracks  of  the  target  ace  patted  in 
Fig.  8,  and  the  corrected  tracks  using  the  fully  developed  bias  errors  are  plotted  in  Fig.  9.  Tho1  results 
show  a  significant  improvement  particularly  in  the  last  half  of  the  flight  but  are  not  as  dramal^tjrs  the 
previous  results  (Fig.  7).  ; 

'*  i 

5.0  SUMMARY 

An  algorithm  has  developed  which  can  be  used  to  detect  and  remove  bias  errors  from  radar 
measurements  when  the  sisasurements  are  made  at  two  separate  and  stationary  radar  sites.  The  algo¬ 
rithm  has  demonstrated  its  ability  to  detect  bias  errors  and  use  this  information  to  converge  tracks  from 
independent  radars.  This  was  accomplished  with  both  real  and  simulated  data.  Since  bias  errors  in 
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Fig.  8  —  Uncorrected  tracks  of  common  target 


Fig.  9  —  Corrected  tracks  of  common  target 

radars  are  the  rule  rather  than  the  exception,  further  investigation  is  essential  if  multiple  platform  sen¬ 
sor  integration  is  to  be  achieved. 
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APPENDIX 


The  position  measurement  vector,  in  a  local  Cartesian  coordinate  system,  can  be  approximated  by 
a  linearized  form  as  follows: 


X 

X 

3  x/Br 

Bx/Br)  Bx/Ba 

A  R  +  N, 

y 

= 

y 

+ 

By  /Br 

By/B-q  By/Ba 

a  tj  +  n. 

Z 

2 

z 

Bz/Br 

Bz/Br)  Bz/Ba 

Aa  +  Na 

where 


(x,y,z)  2 

is  the  target  position  vector  with  respect  to  Site  1  as 
determined  by  measurements  made  at  Site  2 

(x.y.z) 

is  the  mean  or  true  position  vector  of  the  target  with 
respect  to  Site  1 

(Bx/Br,By/Br)  etc.) 

are  the  partial  derivatives  of  the  components  of  the 
position  vector  with  respect  to  changes  in  the  measure¬ 
ments  made  at  Site  2 

( r,r),a ) 

represent  the  range,  elevation,  and  azimuth  measure¬ 
ments  made  at  Site  2 

(A  R,  A7),Aa)2 

represent  bias  errors  in  the  measurements  made  at  Site 
2 

(N„NvNa) 

represents  noise  in  each  of  the  measurements  made  at 
Site  2. 

In  vector-matrix  notation  Eq.  (Al)  can  be  expressed  as 

X2  =  X  +  A2(B2  +  N2).  (A2) 

The  elements  of  the  A  matrix  are  the  partial  derivatives  of  the  components  of  the  position  vector.  The 
components  of  the  position  vector  are  with  respect  to  the  origin  located  at  Site  1,  and  the  partial  deriva¬ 
tives  are  taken  with  respect  to  the  measurement  made  at  Site  2.  The  sites  are  located  in  a 
latitude/longitude  system  in  which  Site  1  has  the  latitude/longitude  pair  of  (Ob,\b)  and  Site  2  has  the 
pair  (0A,\A).  The  parameters  8fl,  8^  represent  the  height  of-  the  respective  sites  above  the  Geoid. 
The  Cartesian  coordinates  of  the  target  with  respect  to  an  origin  located  at  Site  1  are  given  by: 
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The  transformation  from  Site  2  to  Site  1  is  accomplished  with  the  rotation  matrix  Tc  and  the  translation 
matrix  Uc.  They  are  defined  as  follows: 


Tc  = 


cos  (A  ,4  -  XB) 
sin  0g  sin  (A ,4  —  A  5) 

— cos  9  B  sin  (kA  -  Kb) 


-sin  0,4  sin  (A  .4  -  kB) 

sin  0A  sin  0B  cos  (A^  -  kB )  +  cos  0A  cos  9B 

cos  Ga  sin  $B  -  sin  6A  cos  6 B  cos  (A A  —  kB) 


cos  9a  sin  (A a  -kB) 

-cos  0A  sin  9b  cos  (A^  -  kB)  +  sin  9A  cos  9B 
cos0A  cos  9b  cos  (A a  -  Afl)  +  sin  9A  sin  0B 


(A4) 


(a  +  8,4)  cos  0,4  sin  (A^  -  kB) 

-(a  +  8,4)  [cos  0,4  sin  05  cos  (A 4  —  kB)  -  sin  9A  cos; 9B] 

(a  +  8,4)  [cos  9a  cos  0fl  cos  (A ,4  -  kB)  +  sin  9A  sin  0B]  -  (a  +  85) 


The  elements  of  matrix  A2  can  now  be  derived  by  differentiating  Eq.  (A3)  and  recalling  that 

x2  =  r2  cos  172  sin  a 2, 
y2  =  r2  cos  172  cos  a2, 
z2  -  r2  sin  t)2. 


(A5) 


(A6) 


They  are  as  follows: 

dxi/dr2  =  T n  cos  r)2  sin  a2  +  Tu  cos  172  cos  £*2  +  T\3  sin  172 
dxi/d^^  -Tur2  sin  172  sina2-  T\2r2  sin ij2  cos a2  4-  Tl3r2  cost)2 
dx)/da2  =  Tur2  cos  t}2  cos  a2  -  Tl2r2  cos  7)2  sin  a2 
dy\/dr2  =  T2\  cos  r\2  sin  a2  +  T22  cos  172  cos  a2  +  T23  sin  t]2 

dy i/9’J2  =  -  T2lr2  sin  172  sin  a2  -  T22r2  sin  t\2  cos  a2  +  T23r2  cos  7)2  (A7) 

dy i/da2  =  T2lr2  cos  172  cos  a2  -  T22r2  cos  172  sin  a2 

dzj/dr2  =  r31  cos  r)2  sin  a2  +  T32  cos  rj2  cos  a2  +  T33  sin  tj2 

dzi/dt)2=-T3ir2  sinr)2  sin a2  -  T32r2  sin  172  cos (*2  +  T33r2  cos rj2 

dz i/da2  =  7’31r2  cos  t\2  cos  £*2  -  T32r2  cos  rj2  sin  0:2. 

where  Tu  is  an  element  of  the  Tc  matrix.  The  elements  of  the  matrix  A  j  in  Eq.  (1)  are  easily  derived 
from  the  relationships 

xj  *»  rx  cos  17 1  sinai 

^1  ->  r\  cos  171  cosai  (A8) 

zx  =  rj  sin  17, 
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